({OSR-TH-  94  0186' 

DEVELOPMENT  OF  PREDICTIVE  REACTION 
MODELS  OF  SOOT  FORMATION 


AD-A278  366 


Approved  for  public  release  { 
distribution  unlimited* 


Final  Technical  Report 

(Grant  No.  AFOSR  91-0129) 


Covering  the  period  from  January  1,  1991  to  December  31,  1993 


prepared  for 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
Bolling  AFB 


r*  FlCTE 
APR  2 11994 


Michael  Frenklach 
Fuel  Science  Program 

Department  of  Materials  Science  and  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Jhis  document  has  oeen  approve 
for  pubhc  release  and  safofit, 
qistribunoD  is  uniirited 


February  9,  1994 


94-12037 


: 


;TiSD3 


M  4  20  109 


REPORT  DOCUMENTATION  PAGE 


om  m.  ew-ctn 


u  Acuor  uss  mu  INN  j  x  report  date 

•  IFebruarv  9,  1994 


4.  mu  AMO 

(U)  Development  of  Predictive  Reaction  Models  of 
Soot  Formation 


A.  AUTMOKS) 

Michael  Frenklach 


IMPORT  _ 

Final  Technical  Report 


PE  -  61102F 
PR  -  2308 
SA  —  BS 

G  -  AFOSR  91-0129 


ORGANIZATION  NAMEIS)  AND  AOORUKU) 

The  Pennsylvania  State  University 
University  Park,  PA  16802 


MOSR-TR-  9  4  0186 


9.  SPOMSORMG /  MONITORING  ASINCT  NAME(S)  AND  AOOAESSUS) 

AFOSR/ NA 

110  Duncan  Avenue,  Suite  B115 
Bolling  AFB  DC  20332-  0001 


12x  OUTUUTKM/AVAILAItUTY  STATEMENT 

Approved  for  public  release;  distribution  is 
unlimited 


ASINCT  REPORT 


tib.  distribution  coos 


IX  ABSTRACT  (Minmwi dOO  woW) 

The  ultimate  objective  of  the  present  research  program  is  to  develop  a  predictive  reaction  model  for  , 
soot  formation  in  hydrocarbon  flames.  The  principal  accomplishments  in  the  reported  period  are: 

( 1)  A  computer  algorithm  was  developed  that  calculates  optical  properties  of  an  ensemble  of 
particles  whose  size  distribution  is  given  in  terms  of  moments  of  the  size  distribution  function.  (2) 
A  computational  study  of  sooting  limits  in  laminar  premixed  flames  was  initiated  and  completed.  1 1 
was  found  that  the  critical  equivalence  ratios  for  soot  appearance,  both  the  absolute  values  and 
temperature  dependencies,  can  be  predicted  fairly  close  to  the  experimental  observations.  (3)  New 
estimation  techniques  were  developed  and  applied  for  calculations  of  standard-state  enthalpies  of 
formation  and  binary  gaseous  diffusion  coefficients  of  polycyclic  aromatic  hydrocarbons  (PAHs) 
and  their  radicals,  thus  providing  critical  information  for  accurate  modeling  of  soot  formation  in 
flames.  (4)  Theoretical  studies  of  a  bench-mark  ion-molecule  reaction  were  initiated  and  completed. 
The  results  obtained  further  support  the  neutral-species  reaction  pathway  as  the  predominant  route 
for  the  formation  and  growth  of  PAHs,  the  precursors  to  soot  in  hydrocarbon  flames.  (5)  A 
computational  study  of  pressure  effect  on  soot  formation  was  performed.  (6)  A  reduced  model  for 
PAH  and  soot  formation  in  turbulent  reactive  flows  was  developed 
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INTRODUCTION 


The  ultimate  objective  of  the  present  research  program  is  to  develop  a  predictive  reaction 
model  for  soot  formation  in  hydrocarbon  flames.  The  specific  objectives  of  the  proposed  3-year 
study  were:  1)  To  extend  the  modeling  efforts  to  computer  simulation  and  analysis  of  more 
complex  sooting  phenomena,  such  as  sooting  limits  in  laminar  premixed  flames;  and  2)  Further 
refinement  of  the  underlying  reaction  mechanism  of  soot  formation. 

ACCOMPLISHMENTS 

The  principal  accomplishments  of  the  past  three  years,  under  the  support  of  this  grant,  are 
the  following: 

•  A  computer  algorithm  was  developed  that  calculates  optical  properties  of  an  ensemble  of 
particles  whose  size  distribution  is  given  in  terms  of  moments  of  the  size  distribution  function. 
The  new  algorithm  allows  us  to  simulate  light  absorption  and  scattering  without  assuming  a 
functional  form  for  the  particle  size  distribution,  and  thus  to  compare  directly  the  numerical 
predictions  of  the  model  to  the  actual  measured  properties,  such  as  intensity  of  scattered  light 
determined  by  laser  diagnostics. 

•  A  computational  study  of  sooting  limits  in  laminar  premixed  flames  was  initiated  and 
completed.  It  was  found  that  the  critical  equivalence  ratios  for  soot  appearance,  both  the 
absolute  values  and  temperature  dependencies,  can  be  predicted  fairly  close  to  the  experimental 
observations.  Sensitivity  and  reaction-path  analyses  were  performed  to  examine  the  factors 
responsible  for  the  computed  behavior. 

•  New  estimation  techniques  were  developed  and  applied  for  calculations  of  standard-state 
enthalpies  of  formation  and  binary  gaseous  diffusion  coefficients  of  polycyclic  aromatic 
hydrocarbons  (PAHs)  and  their  radicals,  thus  providing  critical  information  for  accurate 
modeling  of  soot  formation  in  flames. 

•  Theoretical  studies  of  a  bench-mark  ion-molecule  reaction  were  initiated  and  completed.  The 
results  obtained  further  support  the  neutral-species  reaction  pathway  as  the  predominant  route 
for  the  formation  and  growth  of  PAHs,  the  precursors  to  soot  in  hydrocarbon  flames. 

•  A  computational  study  of  pressure  effect  on  soot  formation  was  performed. 

•  A  reduced  model  for  PAH  and  soot  formation  in  turbulent  reactive  flows  was  developed. 

The  results  obtained  are  detailed  below;  for  those  that  already  appeared  in  the  published  form 
only  brief  descriptions  are  given,  and  for  those  yet  not  published,  a  more  detailed  account  is 
provided. 

RESULTS 

Detailed  Simulation  of  Sooting  Limits  in  Laminar  Premixed  Flames 1 

The  critical  equivalence  ratios  were  determined  for  several  atmospheric  laminar  premixed 
flames  of  C  2-fuels:  ethane,  ethylene  and  acetylene.  The  value  of  the  critical  equivalence  ratio,  <pc, 
for  a  given  fuel  and  a  given  maximum  flame  temperature,  Tm,  was  determined  by  computing 
flames  with  different  equivalence  ratios.  The  maximum  flame  temperature  in  these  runs  was 
maintained  approximately  the  same  by  adjusting,  similarly  to  the  experiment,  the  amount  of  N2 
in  the  mixture.  For  each  of  the  fuels,  at  least  two  series  of  flames  were  simulated,  each  at  a 
different  maximum  flame  temperature,  in  order  to  test  the  temperature  dependence  predicted  for 
0c 
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The  numerical  simulation  of  soot  particle  formation  was  performed  in  two  stages:  first,  the 
production  of  the  initial  PAH  species;  and  second,  PAH  further  growth  and  nucleation  and 
growth  of  soot  particles.  In  the  first  stage,  PAH  formation  up  to  coronene  was  simulated  in  a 
burner-stabilized  flame  configuration  with  a  specified  temperature  profile  and  flow  rate.  The 
maximum  temperature  reached  in  the  free-flame  simulation,  with  a  given  0,  was  matched  to  the 
corresponding  experimental  value  measured  by  Harris  et  al.2  The  adjustment  of  the  maximum 
flame  temperature  was  accomplished,  similarly  to  the  experiment,  by  varying  the  amount  of  N2 
present  in  the  mixture,  but  not  necessarily  equal  to  the  experimental  level  of  dilution,  and  being 
typically  5%  higher  in  N2  mole  fraction.  In  the  second  stage,  the  profiles  of  H,  H2,  C2H2,  O2, 
OH,  H2O  and  acepyrene  (A4R.5)  obtained  in  the  flame  simulation  were  used  as  input  for  the 
simulation  of  the  growth  of  PAHs  beyond  acepyrene  by  the  technique  of  chemical  lumping  and 
the  soot  particle  nucleation  and  growth  by  a  method  of  moments. 
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Figure  1.  Comparison  of  the  computed  and  experimentally  determined  critical  equivalence  ratios 
0c  for  ethane,  ethylene  and  acetylene  flames,  as  a  function  of  the  maximum  flame 
temperature  Tm.  The  open  symbols  represent  the  experimental  data  of  Harris  et  al.,2 
solid  lines — linear  fits  to  the  experimental  points,  and  solid  symbols — the  present 
computational  predictions. 


The  critical  equivalence  ratios  determined  in  this  manner  for  the  three  fuels  are  displayed  in 
Fig.  1,  where  they  are  compared  with  the  experimental  data  of  Harris  et  al.2  As  can  be  seen  in 
this  figure,  the  computed  critical  equivalence  ratios  follow  the  correct  temperature  dependence  for 
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each  fuel,  i.e.,  the  critical  equivalence  ratio  increases  as  the  flame  temperature  increases,  and  the 
correct  dependence  on  fuel  type,  for  a  constant  flame  temperature  Tm,  i.e.,  acetylene  soots  earlier 
than  ethylene  which  soots  earlier  than  ethane.  The  quantitative  agreement  between  the  computed 
and  experimental  values  of  <pc  is  also  reasonable,  considering  the  lack  of  the  precise  experimental 
data  for  model  input  and  boundary  conditions,  uncertainties  in  the  model  parameters,  and  the  fact 
that  the  reaction  mechanism  was  adopted  from  the  previous  study3  without  any  adjustments. 

The  analysis  of  the  computational  results  strongly  suggested  that  the  appearance  of  soot  and 
hence  the  sooting  limits  in  these  flames  is  controlled  by  two  factors:  concentration  of  acetylene 
and  the  growth  of  PAHs.  The  second  factor — the  PAH  build-up — is  limited  by  the  rise  in  flame 
temperature  towards  the  end  of  the  main  reaction  zone.  The  oxidation  of  PAHs,  the  molecular 
precursors  to  soot  particles,  was  not  found  to  be  a  controlling  factor.  This  does  not  mean, 
however,  that  hydrocarbon  oxidation  does  not  play  a  role  in  soot  formation  On  the  contrary, 
oxidation  of  pre-aromatic-ring  species,  like  C2H3  by  O2,  determines  to  a  great  degree  the 
concentration  level  of  phenyl,  the  first  aromatic  ring.  Also,  the  oxidation  of  PAHs  by  OH  begins 
to  be  more  pronounced  in  locations  away  from  the  main  flame  zone,  where  OH  becomes  the 
dominant  oxidant. 

Evaluation  of  the  Enthalpies  of Formation  of  PAH  Molecules  and  Radicate* 

As  part  of  the  ongoing  efforts  on  soot  formation  model  development  and  refinement,  the 
enthalpies  of  formation  of  polycyclic  aromatic  hydrocarbon  (PAH)  molecules,  radicals  and 
substituted  aromatics  were  investigated.  The  enthalpy  data  used  in  previous  and  current 
modeling  efforts  have  been  mainly  derived  from  Benson’s  group  additivity  method.  Recent 
results  from  ab  initio  quantum  mechanical  calculations  have  shown  that  for  large  peri -condensed 
PAHs,  the  enthalpies  of  formation  predicted  by  the  group  additivity  method  could  be  off  by  as 
much  as  10  kcal/mol.5  The  objective  of  our  work  was  to  determine  accurately  the  enthalpies  of 
formation  for  the  major  aromatic  species  involved  in  the  growth  of  PAHs  and  soot  in  high- 
temperature  environments.  These  species  include  PAH  molecules  up  to  3-circumcoronene,  PAH 
radicals  up  to  coronyl,  and  some  substituted  aromatics. 

A  method  for  accurate  and  economical  estimation  of  the  enthalpies  of  formation  for 
benzenoid  aromatic  species  was  developed  which  combines  semiempirical  quantum-mechanical 
calculations  with  group  corrections.  In  this  method,  the  deviation  between  experimental  and 
calculated  enthalpies  of  formation  is  partitioned  into  structural  groups.  The  general  idea  is  to  use 
a  calculation  method  which  is  computationally  less  demanding  than  quantum  ab  initio  but 
physically  more  realistic  than  group-additivity  calculations,  and  then  to  correct  the  results  using  a 
group-additivity  scheme  for  the  numerical  inaccuracy  of  the  semiempirical  quantum-chemical 
method.  The  advantage  of  this  approach  is  that  it  can  be  applied  consistently  to  both  molecules 
and  radicals.  The  obtained  with  this  method  standard-state  enthalpies  of  formation  compare  well 
with  the  experimental  data  and  with  the  results  of  ab  initio  quantum-mechanical  calculations 
available  in  the  literature.  The  calculated  enthalpies  of  formation  for  large  peri-condensed 
aromatic  molecules  are  extended  smoothly  to  the  limit  of  a  graphite  monolayer. 

The  developed  method  was  also  applied  to  examine  the  stability  of  aryl  radicals  as  functions 
of  molecular  size  and  radical  position,  which  indicates  that  the  strength  of  aryl-H  bond  is 
essentially  independent  of  molecular  size,  but  dependent  on  the  neighboring  geometry  of  the  C-H 
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bonds.  Smaller  aryl-H  bond  strengths  are  predicted  for  the  C-H  bonds  situated  in  the  bay  region 
of  the  aromatic  molecules  and  are  the  results  of  larger  steric  repulsion  between  the  neighboring  H 
atoms. 

Transport  Properties  of  Polycyclic  Aromatic  Hydrocarbon s6 

Diffusion  of  gaseous  PAHs  is  an  important  process  that  significantly  affects  the  results  of 
numerical  simulations  of  soot  formation  in  hydrocarbon  flames.  Although  several  empirical 
correlations  are  available  for  estimation  of  diffusion  coefficients  of  aromatic  compounds,  they  are 
not  easy  to  implement  in  the  computer  codes  now  in  use  by  the  combustion  community.  The 
common  practice  is  to  compute  transport  properties,  including  diffusion  coefficients,  through  the 
parameters  of  Lennard- Jones  12-6  (LJ)  potential,  and  hence  the  LJ  collision  diameters  a  and 
potential  energy  well  depths  e,  specified  for  individual  chemical  compounds,  serve  as  input  to 
flame  codes.  Such  approach  assures  not  only  consistency  among  various  transport  properties, 
such  as  gaseous  viscosities  and  thermal  conductivities,  but  also  consistency  in  evaluation  of 
transport  properties  with  calculations  of  reaction  rates,  such  as  those  encountered  in  coagulation 
of  PAH  species. 

A  method  was  developed  for  systematic  evaluation  of  the  Lennard-Jones  parameters  for 
PAH  compounds,  in  which  correlations  for  these  parameters  are  derived  using  a  group 
contribution  technique  for  critical  temperatures  and  pressures  and  the  Tee-Gotoh- Stewart 
correlations  of  corresponding  states.  The  Lennard-Jones  self-collision  diameters  and  well  depths 
of  29  polycyclic  aromatic  hydrocarbons  were  estimated  using  this  approach  and  are  shown  to 
correlate  with  the  molecular  weights  of  aromatics.  The  correlation  obtained  for  the  self-collision 
diameters  exhibits  the  one-third  power  in  the  molecular  weight,  indicating  that,  due  to  rotation, 
planar  gaseous  PAH  molecules  behave  like  spherical  particles.  The  gaseous  binary  diffusion 
coefficients  of  aromatics  in  common  gases  were  calculated  with  Chapman-Enskog  equation  using 
the  estimated  Lennard-Jones  parameters  and  were  found  to  compare  well  with  the  available 
experimental  data  and  the  predictions  of  one  of  the  most  reliable  empirical  approximations. 

Semiempirical  Quantum-Mechanical  and  RRKM  Studies1  of  Reaction  CgH7+  +  C2H2  —» CioH9+ 

The  title  reaction  is  one  of  the  key  steps  in  the  mechanism  proposed  by  Calcote  and 
coworkers.8  Seven  major  isomers  for  the  CioHg*  adduct  and  1 1  possible  product  channels  were 
identified  in  our  study,  indicating  the  complexity  of  the  reaction.  The  energetics  and  molecular 
parameters  of  the  reactants,  intermediate  species,  products,  and  transition  states  were  evaluated 
using  the  semiempirical  quantum-mechanical  AMI  method.  The  rate  coefficients  and  their 
pressure  and  temperature  dependence  were  then  calculated  using  the  RRKM  theory  with  a  full 
consideration  of  angular  momentum  conservation.  The  microscopic  rate  coefficients  for  the 
reaction  channels  involving  the  dissociation  of  the  energized  complexes  were  determined  with  the 
microcanonical  variational  transition  state  theory.  The  calculations  were  performed  for  a 
temperature  range  of  300-2000  K  and  pressures  of  20  and  760  torr.  These  physical  conditions 
correspond  to  experimental  studies  of  ion  formation  in  flames  .9’10 

The  major  results  of  this  study  can  be  summarized  as  follows.  At  temperatures  below 
1000  K,  the  overall  rate  coefficient  is  equal  to  that  of  the  Langevin  limit,  6x1 014  cm-^moHs*1. 
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However,  at  flame  temperatures  above  1500  K,  where  most  of  the  ions  are  observed  in  laminar 
premixed  flames,  the  rate  coefficients  of  the  reactions  between  PAH  ions  and  acetylene  could  be 
as  much  as  three  orders  of  magnitude  lower  than  the  Langevin  limit.  The  main  reason  for  this  is 
insufficient  collisional  stabilization.  Thus,  it  is  expected  that  many,  if  not  all,  aryl  ion-molecule 
reactions  occurring  in  hydrocarbon  flames  may  have  rate  coefficients  on  the  same  order  of 
magnitude  as  those  of  neutral  aryl  radical-acetylene  reactions.  The  latter  are  typically  on  the 
order  of  1012  cm3mol*1s‘1.  This  result  suggests  that  the  rate  coefficients  used  in  the  reaction 
model  of  ion  formation  proposed  by  O’cote  and  coworkers8  may  be  several  orders  of  magnitude 
too  high. 

Pressure  Effect  on  Soot  Formation 

The  ability  to  predict  soot  formation  at  high  pressures  is  of  considerable  interest  because 
most  practical  combustion  devices  operate  at  elevated  pressures.  Yet,  most  laboratory  flame 
experiments  focusing  on  soot  were  performed  at  subatmospheric  and  atmospheric  conditions. 
Recently,  a  systematic  investigation  of  pressure  effect  on  soot  formation  in  laminar  premixed 
flames  of  ethylene  have  been  undertaken  by  Jander,  Wagner  and  coworkers.11'14  In  their  studies, 
the  pressure  was  varied  from  1  to  100  bar.  The  objective  of  the  present  study  was  to  test 
whether  our  recently  developed  soot  formation  model3’15  is  able  to  account  for  these 
observations.  Specifically,  we  chose  to  simulate  the  10  bar  flames  of  Bonig  et  al. 13  To  determine 
the  effect  of  pressure  on  soot  production,  the  simulation  results  were  compared  to  those  obtained 
for  an  atmospheric  flame. 

The  computational  model  consists  of  the  following  principal  processes:3’15  initial  aromatic 
ring  formation  during  small  hydrocarbon  pyrolysis  and  oxidation,  planar  PAH  growth  through 
the  hydrogen-abstraction-acetylene-addition  (HACA)  reaction  mechanism,  particle  nucleation 
through  coalescence  of  PAHs  into  three-dimensional  clusters,  particle  growth  by  coagulation  and 
surface  reactions  of  the  forming  clusters  and  particles. 

A  recently  updated  500-reaction  and  100-species  mechanism  was  used  to  describe  the  initial 
hydrocarbon  pyrolysis  and  oxidation,  and  PAH  formation  and  growth  up  to  coronene.  The 
formation  of  the  first  aromatic  ring  proceeds  via  the  reactions 


C3H3  +  C3H3 

— > 

benzene 

(1) 

W-C4H3  +  C2H2 

-> 

phenyl 

(2) 

£-C6H4  +  H 

phenyl 

(3) 

W-C6H5 

-> 

phenyl 

(4) 

W-C4H5  +  C2H2 

benzene  +  H 

(5) 

f-C6H6  +  H 

benzene  +  H 

(6) 

W-C6H7 

-> 

benzene  +  H 

(7) 

The  rate  coefficient  of  reaction  (1)  was  assigned  a  value  of  1*10' 1  cm3mol-!s_1,  based  on  an 
analysis  of  benzene  formation  in  subatmospheric  and  atmospheric-pressure  flames.  The  rate 
coefficients  of  reactions  (2)  -  (7)  were  calculated  using  the  RRKM  theory  .16 
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Flame  simulations  were  carried  out  using  the  Sandia  burner  code.17  The  transport  properties 
of  PAHs  were  taken  from  our  recent  evaluation.6  The  computed  profiles  of  H,  H2,  C2H2,  02, 
OH,  H2O  and  the  rate  of  pyrene  formation  were  then  used  as  inputs  for  the  simulation  of  PAH 
growth  beyond  pyrene  and  of  soot  formation.  The  PAH  growth  beyond  pyrene  was  described 
by  a  replicating  HACA  reaction  sequence  given  in  Table  1.  The  rate  coefficients  of  the 
replicating  reactions  listed  in  Table  1  were  estimated  based  on  analogous  reactions  of  one-  to 
three-ring  aromatics.  The  formation  of  soot  particles  was  assumed  to  be  initiated  by  the 
coalescence  of  two  PAH  molecules.  The  subsequent  growth  of  soot  particles  proceeds  through 
particle-particle  and  particie-PAH  coagulation  and  through  surface  reactions.  Particle  coagulation 
was  assumed  to  be  in  the  free-molecular  regime,  with  a  size-independent  enhancement  factor  of 
2. 2. 22  Modeling  of  PAH  growth  beyond  pyrene  and  of  soot  particle  dynamics  was 
accomplished  using  the  method  of  moments. 15 

In  the  present  model,  the  surface  growth  is  assumed3*15  to  occur  due  to  reactions  of  surface 
radicals  with  C2H2.  The  surface  oxidation  is  assumed  to  occur  via  reactions  of  surface  radicals 
with  O2  and  of  the  particle  surface  with  OH.  The  surface  reactions  and  their  associated  rate 
coefficients  are  provided  in  Table  2.  The  rate  coefficients  of  the  heterogeneous  reactions,  S1-S4, 
were  estimated  based  on  analogous  gas-phase  reactions  of  one-ring  aromatics.  The  mechanism  of 
soot  particle  oxidation  by  OH  is  not  well  understood  and  therefore  it  was  described  by  the 
probability  of  reaction  when  an  OH  radical  collides  with  the  particle  surface.  This  probability 
was  taken  to  be  0. 13  from  the  work  of  Neoh  et  al.23 

The  rates  of  the  surface  reactions,  S1-S4,  for  the  zth  particle  were  described  by3*15 


*s  Cg  cc  x§  S,N,t  (8) 

where  k$  is  the  per-site  rate  coefficient,  Cg  is  the  concentration  of  gaseous  species  g,  a  is  the 
fraction  of  surface  sites  available  for  reaction,  whose  value  will  be  discussed  below,  xs  is  the 

number  density  of  surface  sites,  and  S,  and  jV,  are  the  surface  area  and  number  density  of  the  /  th 
particle,  respectively. 

Three  10  bar,  burner-stabilized,  ethylene-air  flames  were  selected  for  the  simulation.  To 
examine  the  influence  of  pressure  on  soot  production,  an  1  bar,  burner-stabilized  flame  was  also 
considered.  The  flame  conditions  are  summarized  in  Table  3.  The  following  experimental  data 
were  supplied  by  Jander:24  temperatures  measured  at  different  distances  from  the  burner  surface 
for  Flames  1-4,  soot  volume  fraction  profiles  in  Flames  2-4,  and  in  Flame  1,  the  mole  fractions 
of  major  flame  products  and  PAH  species  at  a  distance  of  3  cm  from  the  burner  surface. 

The  principal  difficulty  in  the  simulation  of  these  10-bar  flames  is  that  the  temperature 
profile  in  the  main-reaction  zone  is  unknown.  In  all  cases,  measurements  were  made  only  in  the 
post-flame  zone,24  therefore  the  temperature  profile  in  the  main-reaction  zone  had  to  be 
extrapolated  from  the  post-flame  measurements.  For  Flames  1  -  3,  the  distance  between  the 
burner  surface  and  the  first  measured  temperature  point  is  4  to  5  times  the  width  of  the  main- 
reaction  zone  with  the  post-flame  temperature  profile  exhibiting  a  marked  curvature  (Fig.  la). 
These  factors  make  the  reliable  temperature  extrapolation  difficult.  It  was  therefore  necessary  to 
impose  additional  constraints  in  order  to  obtain  the  entire  temperature  profile  required  for  the 
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TABLE  1.  Replicating  reaction  mechanism  of  PAH  growth  beyond  pyrene 


k  = 

=  A  TntE/RT 

No. 

React  iona 

A 

n 

E 

Comments 

0 

:m3mol-1s- 

(kJ  mol-1) 

LI 

Am*  +  C2h2 

— >  A,  +  H 

6.6  (33) 

-5.92 

94.5 

1  atm,  b 

12  (33) 

-5.59 

106.9 

10  atm,  b 

Lla 

A,  +  H 

^  Ay  +  H2 

2.5(14) 

66.9 

c 

Lib 

A,  +  OH 

Ay  +  H2O 

2.1(13) 

19.1 

d 

Lie 

A,*  +  H 

^  A, 

6.4  (20) 

-2.15 

7.9 

1  atm,  e 

2.2  (26) 

-3.67 

24.6 

10  atm,  e 

Lid 

Ay  +  O2 

— >  products 

2.0(12) 

31.3 

f 

L2 

A,*  +  C2H2 

A,  C2H  +  H 

5.1  (38) 

-7.09 

123.7 

1  atm,  b 

2.8  (34) 

-5.81 

132.3 

10  atm,  b 

L2a 

a,c2h+  h 

^  Aj€2H-  +  H2 

2.5(14) 

66.9 

c 

L2b 

a,c2h+  oh 

^  A,C2H*  +  H20 

2.1(13) 

19.1 

d 

L2c 

A,€2H-  +  H 

^  a,c2h 

6.4  (20) 

-2.15 

7.9 

1  atm,  e 

2.2  (26) 

-3.67 

24.6 

10  atm,  e 

L2d 

Aj€2H-  +  02 

— »  products 

2.0(12) 

31.3 

f 

L3 

A,€2H«  +  C2H2 

^  Af+i* 

1.4  (52) 

-11.58 

108.9 

1  atm,  b 

1.1  (44) 

-9.08 

97.9 

10  atm,  b 

L3a 

AI+I  +  H 

^  Aj+p  +  H2 

2.5(14) 

66.9 

c 

L3b 

A,+i  +  OH 

^  A#fi*  +  H20 

2.1(13) 

19.1 

d 

L3c 

Al+1*  +  H 

^  A/+i 

6.4  (20) 

-2.15 

7.9 

1  atm,  e 

2.2  (26) 

-3.67 

24.6 

10  atm,  e 

L3d 

A  i+i*  +  O2 

A,  +  products 

2.0(12) 

31.3 

f 

L4 

A|+i*  +  C2H2 

^  A^2  +  H 

6.6  (33) 

-5.92 

94.5 

1  atm,  b 

1.2  (33) 

-5.59 

106.9 

10  atm,  b 

L4a 

Ah-2+  H 

^  A  ,+?•  +  H2 

2.5(14) 

66.9 

c 

L4b 

A, +2  +  OH 

^  Aj4-2*  +  H2O 

2.1(13) 

19.1 

d 

L4c 

A 1+2*  +  H 

^  A, +2 

6.4  (20) 

-2.15 

7.9 

1  atm,  e 

2.2  (26) 

-3.67 

24.6 

10  atm,  e 

L4d 

Ah-2*  +  O2 

— >  Af+i+  products 

2.0(12) 

31.3 

f 

L5 

A, +2*  +  C2H2 

— >  A  j  +  H 

6.6  (33) 

-5.92 

94.5 

1  atm,  b 

!H=_=== 

1.2  (33) 

-5.59 

10  atm,  b 

a  The  rate  coefficients  for  the  reverse  directions  were  determined  via  equilibrium  constants.  °  Based  on  RRKM 
calculations  for  one-  to  three-ring  aryl  reaction  with  acetylene.16  c  Taken  from  Kiefer  et  al.1*  for  the  analogous 
reaction  of  benzene.  d  Taken  from  Madronich  and  Felder19  for  the  analogous  reaction  of  benzene.  e  Based  on 
RRKM  calculations  for  the  phenyl  and  hydrogen  atom  association  reaction  with  parameters  taken  from  Rao  and 
Skinner.:o  /Taken  from  Lin  and  Lin:i  for  the  analogous  reaction  of  phenyl. 
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TABLE  2.  Reaction  mechanism  of  soot  particle  surface  growth  and  oxidation 


k  =  A  Tnt~E/RT 


No. 

Reaction" 

A 

(cm3s-1) 

n 

E 

(kJ  mol-1) 

SI 

Csoot-H  +  H 

Csoot*  +  H2 

4.2  (-10) 

669 

S-l 

Csoot*  +  H2 

— > 

Csoot-H  +  H 

6.5  (-12) 

390 

S2 

Csoot*  +  H 

-> 

Csoot- H 

1.7  (-10) 

S3 

Csoot*  +  C2H2  -> 

Csoot-H  +  H 

14  (-12) 

0.4 

35.1 

S4 

Csoot*  +  C2 

products 

3.7  (-12) 

31.3 

S5 

Csoot-H  +  OH  — > 

products 

0.13* 

"  Csoot-H  represents  an  aim-chair  site  on  the  soot  particle  surface  and  Csoot*  the  corresponding  radical. 
*  Reaction  probability  upon  OH  collision  with  soot  particle  surface.23 


TABLE  3.  Conditions  of  simulated  burner-stabilized  ethylene-air  flames24 


Flame  no. 

Pressure 

(bar) 

C/O 

Cold  gas  velocity 
(cm/s) 

Maximum  flame 
temperature  (K) 

1 

10 

0.68 

6.0 

1880" 

2 

10 

0.673 

3.0 

1895" 

3 

10 

0.60 

6.0 

2017" 

4 

1 

0.69 

5.9 

1711* 

"  Equal  to  the  adiabatic  flame  temperature.  *  Obtained  by  extrapolating  the  post-flame  temperature  data  and 
assuming  the  maximum  flame  temperature  is  located  near  the  point  where  the  first  soot  volume  fraction 
measurement  is  reported. 
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Mole  fraction  Mole  fraction  Temperature  (K) 


Distance  from  the  burner  (cm) 


Figure  1.  Flame  1.  (a)  lines — assumed  temperature  profiles:  solid  line — the  location  of  Tmax 
is  at  the  onset  of  soot  luminosity,  7^  =  Ta a;  dotted  line — Tmax  =  Tac j  -  70  K.  (b) 
and  (c)  concentrations  of  species  indicated,  computed  with  temperature  profiles 
shown  in  (a);  dashed  line  in  (b) — with  consumption  of  acetylene  by  soot  surface 
reaction.  Symbols  in  (a)  and  (b)  are  experimental  data.24 
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simulations.  It  was  assumed  that  (1)  the  flame  temperature  reaches  maximum  at  the  location 
where  the  soot  luminosity  is  first  seen:  0. 1  cm  above  the  burner  for  Flames  1  and  3  and  0.2  cm 
for  Flame  2,  and  (2)  the  maximum  flame  temperature,  rmax>  is  equal  to  the  adiabatic  flame 
temperature,  7^. 

Unlike  Flames  1  -  3,  temperature  extrapolation  can  be  made  with  a  greater  confidence  for  the 
atmospheric-pressure  flame.  Flame  4.  In  this  flame,  measurements  were  made  relatively  close  to 
the  main-reaction  zone.  In  addition,  the  post-flame  temperature  profile  is  nearly  linear. 
Therefore,  T max  was  determined  by  extrapolating  the  temperature  data  to  the  position  where  the 
first  soot  volume  fraction  measurement  was  reported,  i.e.,  0.3  cm  from  ti.e  burner  surface.  The 
maximum  temperature  of  Flame  4  determined  in  this  manner  is  about  200  K  lower  than  7^ 

An  example  of  the  temperature  extrapolation  for  a  10-bar  flame  is  shown  in  Fig.  la.  The  solid 
lines  in  all  panels  of  Fig.  1  were  obtained  with  the  assumption  that  =  7^  A  sensitivity  test 
was  performed  by  assuming  Tmax  =  Ta<\  -  70  K,  which  is  shown  by  the  dotted  line  in  Fig.  la. 
The  mole  fractions  of  H,  C2H2,  benzene  and  pyrene  computed  using  the  two  temperature 
profiles  are  presented  in  Figs,  lb  and  lc.  The  differences  in  the  H  and  C2H2  concentrations 
calculated  with  the  two  temperature  cases  are  insignificant.  For  benzene  and  pyrene,  the  post¬ 
flame  profiles  are  clearly  affected  by  the  assumption  made  on  7’max.  However,  the  most 
important  properties  for  the  nucleation  of  soot  particles — the  initial  part  of  the  concentration 
profiles  and  their  peak  values — are  essentially  unaffected  by  this  assumption.  The 
computational  features  observed  above  are  the  same  for  all  of  the  flames  considered. 

In  Fig  lb,  a  comparison  of  the  model  prediction  and  the  experimental  data  for  acetylene  in 
Flame  1  is  presented.  It  is  apparent  that  the  experimental  data  are  well  predicted  by  the  model. 
The  agreement  is  particularly  good  in  the  region  where  soot  gains  most  of  its  mass  through 
acetylene  surface  reactions,  i.e.,  between  0.2  to  1  cm  from  the  burner  surface.  The  seeming  over- 
prediction  of  the  acetylene  peak  mole  fraction  is  not  caused  by  the  assumption  made  on  Tmax.  as 
the  maximum  mole  fractions  of  acetylene  computed  using  the  two  temperature  profiles  shown  in 
Fig.  la  are  essentially  the  same.  As  shown  in  Fig.  lb,  the  acetylene  concentration  at  3  cm  from 
the  burner  surface  was  more  closely  predicted  with  the  assumption  Tmax  =  rad  than  using  lower 
values  of  Tmax.  As  will  be  discussed  later,  acetylene  consumption  by  soot  was  found  to  be 
insignificant  over  the  entire  flame 

The  kinetic  model  predicted  reasonably  well  other  measured  species  concentrations,  see 
Fig.  2.  The  computed  main-reaction  zone  width  is  about  0. 1  cm,  which  is  consistent  with  the 
assumed  distance  between  the  burner  surface  and  the  location  of  the  maximum  flame  temperature. 
The  experimental  concentrations  of  H2,  H2O  and  CO2  at  3  cm  from  the  burner  surface  were  fairly 
well  reproduced  by  the  model,  with  the  exception  of  CO  which  was  slightly  over  predicted.  Such 
a  level  of  agreement  is  typical,  and  commonly  seen  for  even  better  studied,  better  spatially 
resolved  sub-atmospheric  flames  (see,  e.g..  Ref.  25). 

As  shown  in  Figs.  2c  and  2d,  the  PAH  concentrations  at  3  cm  from  the  burner  surface  were 
computed  within  a  factor  of  2  to  3  of  the  experimental  values.  For  this  flame,  the  measured  post¬ 
flame  concentrations  of  naphthalene  were  larger  than  those  of  benzene.  This  was  successfully 
predicted  by  the  model.  Analysis  of  reaction  pathways  indicated  that  for  Flames  1  -  3,  benzene 
is  produced  through  reactions  5-7  with  approximately  equal  rates.  With  the  present  assignment 
for  the  rate  coefficient  of  reaction  1,  the  contribution  of  propargyl  recombination  to  benzene 
production  is  negligible. 
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Mole  fraction  Mole  fraction 


Distance  from  the  burner  (cm) 


Figure  2.  Computed  (lines)  and  experimental24  (symbols)  concentrations  of  selected  species 
for  Flame  1. 


A  characteristic  behavior  of  ethylene  flames  is  the  production  of  large  amounts  of  vinyl 
radicals  in  the  preflame  zone.  The  vinyl  radicals  react  quickly  with  ethylene  to  produce  1,3- 
butadiene  (C4H6)  through  C2H4  +  C2H3  -»  C4H6  +  H.  1,3-butadien-l-y!  (/1-C4H5)  is  then 
produced  from  C4H6  through  H-abstraction  by  OH.  The  further  reaction  of  C4H6  with  C2H3 
leads  to  1,3,5-hexatriene  (^-C^g).  1,3,5-hexatrien-l-yl  (w-C^H?)  is  produced  from  ^-C^Hg  by 
the  H-abstraction  reactions.  At  the  same  time,  vinyl  radicals  decompose  to  acetylene  via  the 
reaction  C2H3  — »  C2H2  +  H.  The  addition  of  vinyl  to  acetylene  produces  vinylacetylene  (C4H4), 
which  in  turn  reacts  with  C2H3  to  produce  l,3-hexadien-5-yne  (/-C6H6).  These  reaction 
pathways  are  very  similar  to  those  identified  for  near-sooting-limit  atmospheric  flames  of 
ethylene.1 

The  comparison  of  the  computed  and  experimental  profiles  of  soot  volume  fractions  are 
presented  in  Fig.  3  for  Flames  2  and  3,  and  Fig.  4a  for  Flames  1  and  4.  Inspection  of  these  figures 
indicates  that  the  shape  of  the  experimental  soot  profiles  for  all  of  the  flames  is  well  predicted  by 
the  model.  The  quantitative  agreement  was  obtained  by  adjusting  a  single  parameter — the  fraction 
of  soot  surface  sites  available  for  reaction,  a  in  Eq.  (8). 

Parameter  a  was  introduced3  to  represent  a  steric  factor,  whose  value  should  decrease  with 
the  increase  in  flame  temperature.  This  is  because  at  high  temperatures,  soot  particle  structural 
units  are  expected  to  be  more  mobile  and  align  themselves  in  a  more  concentric  manner,  so  that 
the  access  of  gaseous  species  to  the  reactive  carbon  sites  is  limited.  Figure  5  displays  the  values 
of  a  obtained  in  the  present  and  past15  studies.  It  is  observed  that  a  markedly  decreases  as  the 
maximum  flame  temperature  increases.  Sensitivity  runs  indicated  that  the  numerical  predictions 
for  soot  volume  fraction  are  extremely  sensitive  to  the  assumed  value  of  a.  As  shown  in  Fig.  3a 
for  Flame  2,  a  factor  of  10  difference  in  a  causes  the  soot  volume  fraction  to  differ  by  more  than 
a  factor  of  30.  The  latter  change  in  soot  volume  fraction  is  expected  to  be  larger  than  the 
uncertainties  of  both  the  experiment  and  the  model.  In  addition,  the  absolute  values  obtained  for 
a  lie  in  a  reasonable  range.  Thus,  it  is  likely  that  the  trend  shown  in  Fig.  5  reflects  the 
dependence  of  a  on  flame  temperature.  At  the  same  time,  no  clear  pressure  dependence  of  a  is 
observed  in  Fig.  5. 

An  important  feature  of  the  present  study  is  that  the  experimentally  observed  effect  of 
pressure  on  soot  volume  fraction  is  well  accounted  for  by  the  model  (Fig.  4a).  Experimentally, 
Bonig  et  al.13  determined  in  ethylene-air  flames  that  the  final  soot  volume  fraction, /v«,  is 
proportional  to  P1  for  P  <  10  bar.  Above  10  bar,  the  dependence  of  soot  production  on  pressure 
weakens  and  /v<»  is  proportional  to  P. 14  It  was  suggested14  that  the  weakening  pressure 
dependence  is  due  to  the  increasing  acetylene  consumption  by  soot  surface  reactions,  thus 
depleting  the  principal  surface  growth  species.  While  simulations  over  a  broader  range  of 
pressure  were  not  attempted  in  the  present  study,  the  influence  of  pressure  on  soot  production  is 
examined  by  concentrating  on  the  comparison  between  Flame  1  at  10  bar  and  Flame  4  at  1  bar. 
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Figure  3.  Comparison  of  model  predictions  and  experimental  data24  for  soot  volume  fraction 
in  Flames  2  and  3. 
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Comparison  of  Dames  1  arnM^  »»  i««1>olation 

C/O  equal  «o  0.672  and  0.696. 
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Figure  5.  The  fraction  of  soot  surface  radical  sites  available  for  reaction,  a,  plotted  as  a 
function  of  maximum  flame  temperature,  Tm&x.  The  line  is  a  fit  to  the  data,  a  = 
[tanh(8168/7max  -  4.57)+l]/2. 


Distance  from  the  burner  (cm) 

Figure  6.  Comparison  of  the  acetylene  concentrations  computed  for  Flames  1  and  4. 
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First,  the  degree  of  acetylene  consumption  was  examined  for  Flame  1 .  A  series  of  iterative 
flame  simulations  was  performed.  Each  iteration  included  the  consumption  of  acetylene  by  its 
surface  reaction  with  soot  particle  and  used  the  particle  surface  area  obtained  from  the  previous 
iteration.  The  final,  converged  result  for  the  acetylene  concentration  is  shown  in  Fig.  lb.  A 
comparison  with  the  original  profile  indicates  that  an  observable  amount  of  acetylene  is  indeed 
consumed  by  soot.  However,  the  decrease  in  the  acetylene  concentration  is  quite  small  and  has 
little  effect  on  the  computed  soot  volume  fraction  under  the  conditions  of  Flame  1 .  The  deviation 
between  the  experimentally  determined  post-flame  concentrations  of  acetylene  and  the  computed 
concentrations  is  significantly  smaller  than  that  obtained  by  Bohm  et  al.14 

To  understand  the  factors  governing  the  observed  pressure  dependence  of  soot  production, 
the  rates  of  nucleation,  coagulation  and  surface  growth  of  the  soot  particles  were  examined. 
Inspection  of  Fig.  4b  indicates  that  the  peak  nucleation  rate  of  soot  in  Flame  1  is  about  three 
orders  of  magnitude  higher  than  the  peak  rate  in  Flame  4.  This  is  caused  by  larger  concentrations 
of  PAHs  produced  in  Flame  1.  The  nucleation  is  balanced  by  coagulation  for  both  flames.  The 
resulting  peak  number  density  of  soot  particles  in  Flame  1  is  about  a  factor  of  50  larger  than  the 
peak  value  in  Flame  4.  Due  to  the  large  coagulation  rate,  the  soot  number  density  in  Flame  1 
decreases  quickly,  approaching  the  value  for  Flame  4  after  ~20  ms  of  reaction  time.  During  this 
same  reaction  period,  the  soot  volume  fraction  increases  rapidly  due  to  the  surface  growth 
reactions,  and  nearly  reaches  its  final  value  by  the  end  of  the  reaction  period  Thus,  soot  gains 
most  of  its  mass  during  the  same  period  of  intense  coagulation. 

The  rates  of  surface  growth  differ  in  Flames  1  and  4  by  approximately  two  orders  of 
magnitude  during  the  initial  reaction  times  when  soot  gains  most  of  its  mass.  Since  surface 
growth  provides  the  ultimate  soot  mass,  the  quantitative  difference  in  between  Flames  1  and  4 
is  caused  by  the  102  difference  in  the  surface  growth  rates.  The  soot  surface  growth  under  the 
conditions  of  the  simulated  flames  is  predominantly  due  to  the  reactions  of  acetylene  with  the 
soot  particle  surface.  The  surface  growth  rate  is  thus  proportional  to  the  acetylene 
concentration.  In  the  region  where  soot  gains  most  of  its  mass,  the  absolute  concentrations  of 
acetylene  in  Flames  1  and  4  differ,  due  to  the  pressure  difference,  by  only  one  order  of  magnitude 
(Fig.  6).  Thus,  the  difference  in  the  acetylene  concentration  alone  cannot  fully  explain  why  the 
surface  growth  rate  in  Flame  1  is  two  orders  of  magnitude  higher  than  that  in  Flame  4. 

Soot  surface  area  is  an  additional  factor  contributing  to  the  growth  rate.  The  surface  area  is 
determined  not  only  by  surface  growth,  but  also  by  nucleation  and  coagulation.  An  additional 
test  was  performed  to  determine  the  contribution  of  nucleation  to  /v„ ».  In  this  test,  the  peak 
nucleation  rate  of  Flame  1  was  artificially  reduced  to  1014  cm-3s-1,  i.e.,  equal  to  that  of  Flame  4. 
The  resulting  was  lowered  by  about  an  order  of  magnitude.  Thus,  the  larger  nucleation  rate 
and  the  higher  acetylene  concentration  contribute  about  equally  to  the  increased  value  of /v»  in 
Flame  1.  The  above  analysis  suggests  that  for  soot  formation  under  the  conditions  simulated 
nucleation  is  as  critical  as  surface  growth. 
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As  described  earlier,  the  present  model  assumes  that  the  nucleation  of  soot  particles  is 
initiated  by  PAH  coalescence.  Consequently,  the  nucleation  rate  depends  strongly  on  PAH 
concentrations.  Inspection  of  the  computed  results  indicates  that  the  influence  of  pressure  on  the 
total  concentration,  particularly  on  the  concentration  of  acetylene,  has  the  strongest  effect  on  the 
production  of  PAHs  and  hence  on  the  nucleation  of  soot  particles. 

A  secondary  effect  of  pressure  on  the  production  of  PAHs  arises  from  the  pressure 
dependence  of  reaction  rate  coefficients.  This  pressure  effect  arises  from  the  reactions  whose 
energized  adducts  undergo  collisional  stabilization.  A  test  run  under  the  conditions  of  Flame  1, 
but  with  the  pressure-dependent  rate  coefficients  set  equal  to  those  at  1  bar,  results  in  the  peak 
concentration  of  naphthalene  that  is  a  factor  of  4  lower  than  that  computed  with  the  rate 
coefficients  at  10  bar. 

It  is  thus  demonstrated  that  the  detailed  chemical  kinetic  model  of  soot  formation,  developed 
and  tested  earlier  for  subatmospheric  and  atmospheric  flames,  can  quantitatively  reproduce  the 
experimental  soot  volume  fraction  profiles  collected  in  a  number  of  10  bar,  burner-stabilized, 
ethylene-air  flames.  Analysis  of  the  computational  results  suggests  that  the  consumption  of 
acetylene  by  soot  particle  surface  growth  is  insignificant  up  to  10  bar.  In  addition,  the  influence 
of  pressure  on  soot  production  arises  mainly  from  a  larger  concentration  of  acetylene  at  high 
pressure,  which  affects  both  the  soot  surface  growth  and  the  nucleation  of  soot  particles. 

Reduced  Model  for  PAH  and  Soot  Formation 

The  entire  reaction  mechanism — from  the  initial  hydrocarbon  to  soot  particles — can  be 
divided  into  three  parts,  each  reduced  with  a  different  method:  small-molecule  reactions 
describing  the  main  combustion  environment  and  the  formation  of  the  first  aromatic  ring;  its 
growth  to  pyrene;  and  nucleation,  growth  and  oxidation  of  soot  particles.  The  former  part,  small- 
molecule  chemistry,  is  a  subject  of  continuous  efforts  of  the  research  community;  and  the  latter, 
particle  dynamics,  has  been  addressed  by  us  previously  using  the  method  of  moments.  Here,  a 
reduced  model  for  PAH  formation  and  growth,  the  central  part  of  the  soot-precursor  chemistry, 
is  presented. 

A  reaction  sequence  leading  from  benzene  to  pyrene,  adopted  for  the  present  study,  is  shown 
in  Fig.  7.  This  reaction  scheme  is  based  on  both  our  past  results1-3’15  and  a  reaction  flux  analysis 
performed  in  the  present  study  using  an  extensive  reaction  set  under  the  conditions  typical  of 
high-pressure  turbulent  combustion.  The  nomenclature  for  chemical  species  used  in  Fig.  7  and 
throughout  the  manuscript  has  been  introduced  previously.26  Briefly,  A/  represents  a  species 
composed  of  /  fused  aromatic  rings;  e  g.,  Ai  is  benzene  and  A4  pyrene.  A/C2H  is  an  A,  aromatic 
with  one  of  its  H  atoms  substituted  with  C2H.  A,--  denotes  a  radical  formed  from  A,  and 
A/C2H*  a  radical  of  A/C2H  with  the  unpaired  electron  located  next  to  the  C2H  group.  P2 
represents  biphenyl,  an  aromatic  molecule  with  non-fused  rings,  whose  importance  was  predicted 
for  the  conditions  rich  in  benzene.27  As  will  be  demonstrated  later  in  the  text,  the  P2-channel 
plays  a  key  role  under  the  conditions  tested  in  the  present  study. 
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Figure  7.  Reaction  scheme  of  PAH  growth. 
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TABLE  4.  Mechanism  of  PAH  growth. 


No. 

Reaction" 

k  = 

A 

(cm3mol~1s-1) 

Ar*™7 

n  E 

(kJ  mol-1) 

References/ 

Comments 

79/ 

A,-  +H  ^  A,-  +  H2 

2.50(+14) 

66.9 

ej 

80/ 

A,  +  OH  ^  At  +  H20 

2.10(+13) 

19.1 

e,g 

81 1 

A,-  +  H  A, 

5.80(-12) 

8.29 

-33.8 

h,i 

82/ 

A,  +  OH  — »  products 

1.00  (+11) 

83/ 

A,-  +  O2  — »  products 

2.00  (+12) 

31.3 

h,  k 

84 

Ai-  — » linear  products 

1.04  (+29) 

-4.26 

324.0 

l 

85 

Ar  +  C2H2  ^  AiC2H  +  H 

2.34  (+34) 

-5.71 

129.4 

m 

86 

A2-  +  C2H2  -  A2C2H  +  H 

4.63  (+31) 

-5.05 

123.0 

m 

87/ 

A,02H*  +  C2H2  ^  A/+i- 

3.31  (+50) 

-10.95 

112.8 

m,  n 

88 

A3-  +  C2H2  — >  A4  +  H 

1.27  (+44) 

-8.71 

138.2 

m 

89 

Ai  +  Ai-  P2  +  H 

1.30  (+03) 

2.72 

11.3 

0 

90 

P2'  +  C2H2  ^  A3  +  H 

1.27  (+44) 

-8.71 

138.2 

^90=^88 

a  The  rate  coefficients  for  the  reverse  directions  were  determined  via  equilibrium  constants. 
e  A,  represents  Ai,  A1C2H,  A2,  A2C2H,  A3  and  P226  f  Taken  from  Kiefer  et  al.18  for  the 
analogous  reaction  of  benzene.  g  Taken  from  Madronich  and  Felder19  for  the  analogous 
reaction  of  benzene.  h  A,-  represents  A]-,  A1C2H*,  A2-,  A2C2H*,  A3-  and  P2-.26  '  Based 
on  RRKM  calculations  for  the  phenyl  and  hydrogen  atom  association  reactions  at  10  bar 
with  the  parameters  taken  from  Rao  and  Skinner.20  The  forward  rate  coefficient  k  = 
A 7" exp (-E/RT -  2.77xl0-3  T -  5.58xlO-7  P-).  J  For  /  =  1  a  linear  hydrocarbon  product  is 
assumed;  for  1  >  1  the  product  is  an  aromatic  molecule  with  two  carbon  atoms  less  than  A,. 
The  rate  coefficient  is  estimated,  see  Ref.  1 .  k  For  i-  1  a  linear  product  is  assumed;  for  i  >  1 
the  product  is  an  aromatic  radical  with  two  carbon  atoms  less  than  A,.  The  rate  parameters 
were  taken  from  Lin  and  Lin21  for  the  analogous  reaction  of  phenyl.  1  Assigned  the  rate 
coefficient  of  the  ring  opening  reaction  Aj-  — » /1-C6H5.  m  Based  on  RRKM  calculations  at 
10  bar  with  parameters  taken  from  Ref.  16.  ni  =  1,  2.  0  The  rate  coefficient  expression  was 
obtained  using  RRKM  calculations,  fitted  to  the  experimental  data  of  Fahr  and  Stein28  and 
Scaiano  and  Stewart.29 
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Assuming  steady  state  for  all  the  radicals  in  the  sequence  shown  in  Fig.  7  and  partial 
equilibrium  for  reactions  (85)  and  (86)  (see  Table  4),  we  obtain  an  analytical  expression  relating 
the  concentrations  of  pyrene  and  benzene, 

““  =  A88[C2H2]^[AlinBi  +  B,E  +  d) 


where 

A  =  MC2H2]{*_87i2  +  *88[C2H2]  +  *83[  02]  +  Z25it82[OH)/(z1>5  +  &82[OH])}  1 
&2  =  ^1,4  (^2,4  +^87[C2H2]  +  ^83[02]-/:_87,251) 


%  =  k86[C2H2]/{k_86[H}) 


^4  “  ^87  [C2H2  ](^_87,  j  +A:85[02]  +  Z2  3) 


^5  “  +*87[C2H2]  +  ^3[02]-*_87il54) 

%  =  *85[C2H2]/(fc_85[H]) 

Bj  =  ZU{Z2A  +  k82  +  k83[02]  +  Zl6  k89[Ax]/[zl6  +  fc_89[  H]  +  Z2>6  fc_89[H]/(/:88[C2H2l)]}-1 
C  =  Zl>6B7*:88*89[C2H2][A1]2{z1,6fc88[C2H2]+  fc_89[H](Z2,6  +  ^8[C2H2  ])}~£  +  A:82[OH][A4] 
D  =  Zl5c{(fc_87>2  +  *88[C2H21  +  fc83[02])(z15  +  fe82[OH])  +  Z2'5k82[OH]Y' 

E  =  ^_87,2^{^2,4  +  ^87CC2H21  +  ^83[02]  -  &_87,25i} 

Zic  =  ^79,j[H]  +  ^8o,i[OH]  +  /:_8[, 


Z2.i  =*-79,,[H2]  +  *-80|J[H2O]  +  *8i,i[H]  . 


The  rate  coefficients  appearing  in  the  above  equations  are  given  in  Table  4  (reactions  79  -  90). 
They  were  obtained  from  RRKM  calculations  for  a  pressure  of  10  bar  and  temperatures  from 
500  to  3000  K.  The  rate  coefficients  evaluated  at  10  bar  are  within  a  factor  of  2  from  the  values 
at  100  bar  and  were  therefore  assumed  to  be  pressure  independent  in  the  present  calculations. 


23 


The  accuracy  of  the  reduced  model  of  PAH  growth  was  tested  by  performing  kinetic 
calculations  on  an  individual  fluid  cell  whose  initial  conditions  are  typical  of  those  producing  soot 
in  high-pressure  turbulent  combustion.  The  results,  shown  in  Fig  8,  indicate  that  the  reduced 
model  overpredicts  the  pyrene  concentration  in  the  beginning,  but  gets  closer  to  it  at  times  around 
1  ms,  the  average  cell  reaction  time.  The  difficulty  of  obtaining  a  more  accurate  reduced  model  for 
the  reaction  scheme  in  Fig.  7  resides  in  the  complex  behavior  of  A2  and  A3.  Assuming  their 
formation  irreversible  and  omitting  the  P2  channel  predicts  well  the  initial  part  of  the  A4  profile 
but  greatly  underpredicts  it  in  the  range  of  interest.  Assuming  A2  is  in  a  steady  state  or 
neglecting  A2  altogether  overpredicts  the  A4  concentration  at  all  reaction  times.  The  physical 
reason  for  this  is  the  compound  kinetics  of  A2:  it  accumulates  at  first  and  then  reacts  at  the  later 
stages  of  reaction.  A  compromise  solution  was  obtained  by  assuming  the  formation  of  A2 
irreversible,  neglecting  the  formation  of  A1C2H  in  the  oxidation  of  A2  by  OH,  and  including  the 
P2  channel,  whose  importance  in  this  case  was  evident  (compare,  e  g.,  curve  b  with  a  in  Fig.  8). 
In  doing  so,  species  P2,  P2-  and  A3  were  assumed  to  be  in  steady  state. 


Figure  8.  Calculated  mole  fraction  of  A4  for  the  conditions  of  a  single  fluid  cell:  Po  =  39  atm, 
To  =  1954  K,  0=  3.68,  and  initial  mole  fraction  of  A4  9.3xl0’12.  a)  Computed 
with  the  full  reaction  scheme  for  PAH  growth  (Fig.  7  and  reactions  79  -  90  in 
Table  4).  b)  The  same  as  (a)  but  without  the  P2  reaction  channel  (i.e.,  without 
reactions  89  and  90).  c)  Computed  with  the  reduced  model  of  PAH  growth. 
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PRINCIPAL  CONCLUSIONS 


The  principal  accomplishment  of  the  last  three-year  period  is  the  demonstration  that  the 
kinetic  model  developed  for  soot  formation  is  capable  of  quantitative-level  predictions  for  not 
only  low-pressure  laboratory  flames,  but  also  for  more  “practical”  flames:  sooting  limits  of 
laminar  premixed  flames,1  soot  profiles  in  high  pressure  flames,30  and  soot  formation  in  turbulent 
combustion.31 

New  theoretical  results  provide  further  evidence  against  an  ionic  pathway  for  polycyclic 
aromatic  precursors  to  soot  in  hydrocarbon  flames. 

New  methods  for  estimation  of  thermodynamic  and  transport  properties  of  polycyclic 
aromatic  hydrocarbons  were  developed. 

FUTURE  WORK  DIRECTION 

Based  on  the  success  in  modeling  soot  formation  at  high  pressures,  the  direction  for  the  future 
research  under  the  currently  renewed  AFOSR-sponsored  project  is  to  attempt  modeling  of  soot 
formation  in  environments  even  closer  related  to  practical  combustors,  like  diffusion  and 
turbulent  flames,  along  with  further  refinement  of  the  underlying  reaction  mechanism  of  soot 
formation. 
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